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A Simple Algorithm for Computing Canonical Forms 
H. Ford, L.R. Hunt, Renjeng Su 

Abstract 

It is well known that all linear time-invariant controllable 
systems can be transformed to Brunovsky canonical form by a 
transformation consisting only of coordinate changes and linear 
feedback. However, the actual procedures for doing this have 
tended to be overly complex. The technique introduced here is 
envisioned as an on-line procedure and is inspired by George Mayer's 
tangent model for nonlinear systems. The process utilizes Meyer's 
block triangular form as an intermediate step in going to 
Brunovsky form. The method also involves orthogonal matrices, 
thus eliminating the need for the computation of matrix inverses. 

In addition, the Kronecker indices can be computed as a 
by-product of this transformation so it is not necessary to 

know them in advance. 
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A Simple Algorithm for Computing Canonical Forms 
H. Ford*, L.R. Hunt*, and Renjeng Su** 


I. Introduction 


In his work at NASA Ames Research Center, George Meyer is 
applying the theory of transformations of nonlinear systems to 
linear systems in order to design automatic flight controllers 
for vertical and short take off aircraft [1], [2], [3], [4], [5], [6]. 
In these articles he introduces a linear system (called the «. 
tangent model) which approximates the nonlinear system, and he 
mentions the importance of taking the tangent model to Brunovsky 
[7] canonical form. Thus we feel that an on-line procedure for 
transforming a controllable linear system to Brunovsky form is 
highly desirable. Additionally in [8] the authors together with 
George Meyer present a modification of the tangent model in which 
the procedure introduced here plays a central role in the linearization 


♦Research supported by NASA Ames Research Center under grant NAG2-189 
and the Joint Services Electronics Program under ONR contract 
N0014-76-C1136. 

♦♦Research supported by NASA Ames Research Center under grant 
NAG2-203 and the Joint Services Electronics Program under ONR 
contract N0014-76-C1136 . 
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process itself. This process involves a Taylor Series expansion 
using lie derivatives. 

Our computation of the transformation proceeds in two 

steps: 

1) An orthogonal coordinate change is used to move the linear 
system to Meyer's block triangular form (see [5]). 

2) Once we are in block triangular form, the process becomes 
trivial and formally involves "Lie differentiation" of 
certain coordinate functions. The on-line procedure for 
doing this and its application to automatic flight control 
are given in [5] . 

In step 1) it is not necessary to calculate a matrix inverse and no 
systems of linear equations need be solved. As a by-product of 
step 1) we find the Kronecker indices of the system. 

Understanding the fact that the block triangular form is a 
natural intermediate step in transforming to Brunovsky form is 
our main contribution in this paper. We have computer programs to 
carry out the entire process. 

In section II, we describe how to transform a single control 
linear system to a string of integrators. In section III, we show 
how to generalize the results of section II to a multi-control 
system. Detailed algorithms for transforming linear systems to 
Brunovsky form are given in section IV, including pseudocode 
programs. Results achieved using a three control system and 
concluding comments are given in section V. We want to emphasize 
the simplicity and ease of implementation of the algorithm. 
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Though the explanation of why it works may seem cumbersome, the 
algorithm is very straightforward. 


II. Single Control Case 


Definition 2.1 . An n dimensional single control system x = f(x,u) 
is called block triangular if x^ is a function only of x^ , , . . . , x^ + ^ 
for i = 1, ...,n where x n+ ^ = u (see [5]). 

For a linear single control system x = Ax + bu, this is 
equivalent to saying that the square matrix H = [oS] is a lower 
Hessenberg matrix [9], that is, all elements above the first super 
diagonal are zero. Notice that all elements on the first super 
diagonal must be nonzero if we are to have a controllable system. 

The outstanding characteristic of the block triangular system is 
that if we start with x^, it is necessary to take n derivatives 
of x^ with respect to t before reaching the control. 


Definition 2.2. An n dimensional single control system x = f(x,u) 
is called a string of integrators if x^ = x^ + ^ for i = l,...,n 


where x . . = u. 
n+1 


Once we have a block triangular linear system x = Ax + bu, we can 
transform it to a string of integrators y = Ay + £v by simply 
letting y 1 = x x , Y 2 = Y]/ • • • >Y n ~ and v = Y n - 


Theorem 2.1. The above transformation from block triangular form 
to a string of integrators consists of only coordinate changes 
and linear feedback. 


Proof: We know that y^ = x^ as mentioned above. 



we have 


Assuming y. = E a. x. for some i < n - 1 

3=1 3 3 


i i j+1 

= = E a . x . = E a . E a_. 

j=l 3 3 j = 


i+i -j -j “j f k ~k 


i+l 

x,. = E 6, x, 

H=1 


k k 


Thus the state variables y are functions only of the state variables 

n 

x. The new control v = y = E 8, x, 

n k=i k k 


n-1 


n-1 


k+1 

= E 8. + 8_ x = E 8 E a v .x. + 8 [ E a . x. + bu 

k=1 k k n n k , x k 0=1 k,j j p n^ =1 n, 3 j 

n 

= E y x + Y +1 u, 

1=1 l i n 1 

and our feedback is linear. [^] 

If we start with a system x = Ax + bu which is not block 
triangular, we first form the augmented matrix H = |p q] • From 

Stewart [9] it is known that there exists an orthogonal 


n 


transformation C 


CHC' 


C 0 

IP 1 


such that CHC is lower Hessenberg. 


fr 7} r c\ 

”~T 

— 


~ ~T - 1 

CO W 

C 

0 


CAC 

Cb 

Ip ij [o o| 

0 

1 


0 

0 


Thus the transformation y = Cx yields a block triangular system 
y = CAC$r + Cbu. 


Notice 
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C 0 


X 

_ 

Cx 

L u J 

L U J 1 

0 1 


U 


u 


is a transformation 


involving coordinate changes only. Once we have a block triangular 
system, we go to a string of integrators as before. 


III. Multi-Control Case 


For a multi-control system, the situation is slightly more 
complicated than for the single control case. Below is the usual 
definition of Kronecker indices [7] . 
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Definition 3.1 . Let x = Ax + Bu be a time-invariant n dimensional 

linear control system with m controls. 

Let r = rank B 
o 

Tj = rank {B,AB, ... , B}-rank{ B,AB, . . . ,A^ ^B} 

We define the Kronecker indices k. . as the number of r.'s that are 

1 m ^ 

> i. Notice k. > k_ > ... > k and .Z.K. = n. 

- 1-2- -m ]=1 ] 

Definition 3.2. By the Brunovsky [7] canonical form , we mean a 

, . . A A A 

linear system y = Ay + Bv such that A equals 


m 


0 10 . 

0 0 10 


0 0 
0 0 


0 0 


0 0 


0 0 


0 

0 


0 0 
0 0 


1 

0 


0 0 


0 I 0 0 


0 I 0 0 


• \ 


0 10 0 


0 

0 


0 10 .. 

0 0 10 . 


0 

0 


1 

0 


0 0 


0 0 0 


0 0 


0 0 


0 10 . 
0 0 10 


, 00 .. 


. . 0 
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That is, the Brunovsky form consists of m strings of integrators 
whose respective lengths are the Kronecker indices. For instance 
a ten dimensional system with three controls and Kronecker indices 
5,3, and 2 will have Brunovsky form 



y 

y 

y 

y 

y 

y 

y 

y 

y 

y 



ooo 
0 0 0 
0 0 0 
0 0 0 
10 0 
0 0 0 
0 0 0 
0 10 
0 0 0 
0 0 1 





Notice that for a system with only one control, the Brunovsky form 
is a string of integrators. 

Though not necessary for the implementation of our algorithm, 
we shall for theoretical reasons want to be able to identify 
individual Kronecker indices with individual controls. For that 


reason we shall introduce the following alternate definition of 
Kronecker indices. 


Definition 3.3. Let x = Ax + Bu be an n dimensional linear system 

with m controls and with B = [b, , ...,b ]. Consider the vectors 

2 2 

b ,b , , . . . ,b. , Ab , . . . ,Ab, ,A b , . . . ,A b, , - . . until we come to a vector 
m m-1 1 m 1 m 1 

K • 

dependent on the previous vectors, call it A 1 b^. Then is the 



Kronecker index associated with the control u ^. If we continue 
in this manner we will get m Kronecker indices whose sum is n. 

This definition will give us the same Kronecker indices as before. 
However, the subscripts will now associate each Kronecker index 
with a control rather than ordering the Kronecker indices (we can 


obviously renumber our controls so that k, > k_ ><_>...> k ) . 

x— z— j— — m 

This definition will be very useful in showing that the block 
triangular form of a system contains all information on the 
Kronecker indices. 

What do we mean by the block triangular form of a linear 
system with m controls ? The most obvious choice would be a 
system consisting of m block triangular systems, each in one 
control. For instance, a ten dimensional block triangular system 
with three controls and Kronecker indices 5,3, and 2 might be of 
the form 


(3.2) 


V 


X 

X 

0 

0 

0 

0 

0 

x 2 


X 

X 

X 

0 

0 

0 

0 

*3 


X 

X 

X 

X 

0 

0 

0 

*4 


X 

X 

X 

X 

X 

0 

0 

*5 

- 

X 

X 

X 

X 

X 

0 

0 

X 6 


0 

0 

0 

0 

0 

X 

X 

X 7 


o' 

0 

0 

0 

0 

X 

X 

x 8 


0 

0 

0 

0 

0 

X 

X 

X 9 


0 

0 

0 

0 

0 

0 

0 

*10 


0 

0 

0 

0 

0 

0 

0 


0 0 0 


X 1 


0 0 0 

0 0 0 


X 2 


0 0 0 

0 0 0 


x 3 


0 0 0 

0 0 0 


X 4 


0 0 0 

0 0 0 


X, 


X 0 0 



5 

4- 


0 0 0 


X 6 


0 0 0 

x o 0 


X 7 


0 0 0 

X 0 0 


X 8 


0X0 

0 X X 


X 9 


0 0 0 

0 X X 


1 X, „ 


0 0 x 

— J 


L J- Uf ( 

L- J 



Here each of the 0's represents zero. Each of the rightmost X's 


represents a nonzero element. The other X's may or may not 
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be zero. We emphasize that 5 derivatives of x^ must be taken to 
reach a control, 3 derivatives of Xg, and 2 derivatives of x^ , 
precisely the Kronecker indices. 

We found it convenient to rename the state variables, 
ordering them in terms of their distance from the controls. By 
distance, we mean the number of derivatives of a variable we must 



We must take 5 derivatives of x^ to reach the controls, 3 

derivatives of x^ , and 2 derivatives of x^. Remarkably, the only 

thing which affects the number of derivatives necessary to reach 

the controls is the rightmost nonzero element in each row. 

For instance, x„ is a function of x c but not a function of x. for 

** b 3 

j > 6 n<5r u. for 1< i< 3 . is a function of x n but not a function 
of x^q nor u^ for l<_i<_3.. This implies x^ is a function of x^ 
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but. not a function of x^g nor u^ for 1 <_i <_3. is a function 

of U£ but not u^. This implies x’^ is a function of but not u^ 
Notice in determining that it takes three derivatives of x^ 
before we reach the controls, all we used was knowledge of the 
rightmost nonzero elements. Thus a block triangular system with 
ten dimensions, three controls, and Kronecker indices 5,3, and 2 
could have the form 


V 


"x 

X 

0 

0 

0 

0 

0 

0 

0 

0 


“ x l" 


"0 

0 

0 

X 2 


X 

X 

X 

0 

0 

0 

0 

0 

0 

0 


x 2 


0 

0 

0 

X 3 


X 

X 

X 

X 

X 

0 

0 

0 

0 

0 


X 3 


0 

0 

0 

x 4 


X 

X 

X 

X 

X 

X 

0 

0 

0 

0 


x 4 


0 

0 

0 

x 5 

= 

X 

X 

X 

X 

X 

X 

X 

X 

0 

0 


x 5 


0 

0 

0 

x 6 


X 

X 

X 

X 

X 

X 

X 

X 

X 

0 


x 6 


0 

0 

0 

x 7 


X 

X 

X 

X 

X 

X 

X 

X 

X 

X 


x 7 


0 

0 

0 

00 


X 

X 

X 

X 

X 

X 

X 

X 

X 

X 


x 8 


X 

0 

0 

x 9 


X 

X 

X 

X 

X 

X 

X 

X 

X 

X 


X 9 


X 

X 

0 

*10, 


i 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 


x 10 


X 

_ 

X 

X 



In each row, the X furthest to the right is nonzero. The other X's 
may or may not be zero. Notice it still takes 5 derivatives of x^ 
to reach the controls, 3 derivatives of x^ , and 2 derivatives of x^. 

The above discussion motivates our definition of block 
triangular for a linear system with several controls. 


Definition 3.4. Let x = Ax + Bu be an n dimensional system with m 
controls. We say the system is block triangular if the square 
matrix H = ft SI is a generalization of a lower Hessenberg matrix that 



we simply call generalized lower Hessenberg and which we now define. 
First, all elements of H above the mth super diagonal are zero. 
Second, if H(I,J) , with I< n, is a zero element with no nonzero 
element to its right in row I, then all elements above H(I,J) in the 
same super diagonal must also be zero. For instance, 

'xxooooooooooo 

XXXOOOOOOOOOO 
xxxxxoooooooo 
xxxxxxooooooo 
xxxxxxxxooooo 

XXXXXXXXXOOOO 

xxxxxxxxxxooo 
xxxxxxxxxxxoo 
xxxxxxxxxxxxo 
xxxxxxxxxxxxx 
0000000000000 
0000000000000 
0000000000000 


(3.5.) H = 


where the rightmost X's are nonzero, represents a ten dimensional 
block triangular system with three controls and Kronecker indices 
5,3, and 2. Notice that since H(4,7) = 0, we have 
H(3,6) = H(2,5) = H(l,4) = 0. Similarly, H(2,4) = 0 implies 
H (1 , 3) = 0. 

. If x = Ax + Bu is an n dimensional system with m controls which 
is not in block triangular form, we can easily make it block 
triangular in the following way. Form the augmented matrix H = ^ J , 


r c oi 

and find an orthogonal matrix C = Q 1 (where I is the m dimensional 


identity) such that CHC is generalized lower Hessenberg. We will 


explain in the next section precisely how to find C. This matrix 

-~T 

m* ■ r ’ i t i i zi u i i 

C satisfies 


chcT - rs ?] [o a r ?] = [t t c s] 


Thus the transformation y = Cx yields a block triangular system 
y = £ac y + CBu. As in the single control case, this transformation 
is just a change of coordinates on our state space. 


Definition 3.5. By zero pattern we shall mean the pattern of the 
rightmost nonzero elements of the matrix H = qJ of a system 
x = Ax + Bu in block triangular form. That is, for each block 
triangular form there is an n tuple of integers 

so that is the column number of rightmost nonzero element in 

. , • th 

the 1 row. 

I 

Theorem 3.1 . Let x = Ax + Bu be an n dimensional system with m 
controls. There is a onO to one correspondence between the 
possible ordered sets of Kronecker indices and the possible zero 
patterns of the block triangular form. That is, a given set of 
Kronecker indices with given associations with the controls (see 
def. 3.3) results in a distinct and unique zero pattern. 


Comment: This implies we can retrieve the Kronecker indices of a 

system knowing only the zero pattern of the block triangular form. 


Proof : 


Consider the original system x = Ax + Bu. 

[S ?][o o] [q T S] = [o o] be the block triangular represent- 


Let H 
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ation. Notice that y = Ay + Bu has the same Kronecker indices as the 
original system x = Ax + Bu. To see this consider that 

. np _ 

] . 


[ / • ■ • f b^ ; Ah> m , « . . , Ab ^ / • • • ] — [ Cb^ , « • . , Cb , CAC Cb^ , . . • , C AC Cb ^ , 


£ Cb , . . . ,Cb. , CAb , . . • , CAb — C[b , . . . , b . , Ab , * . . / Ab . / . 
m 1 m ± m 1 m 1 

Since C is nonsingular, the rank of any selection of columns from 

[b , ...,b.,Ab , . . . , Ab. , . . . ] is the same as the rank of the 
m 1 m 1 

corresponding selection of columns from [b , . . . ,b^,Ab m , . . . , Ab^ , . . . ] . 
Therefore, we can assume without loss of generality that the 
original system is already in block triangular form. 

The rest of the proof depends on two basic principles from 
linear algebra: 

1. When we multiply a matrix times a vector, the product 
is a linear combination of the columns of the matrix. 

2. If we have a collection of n nonzero vectors v. , . . . ,v , 

1 n 

then the linear combination with a^O, is 

linearly independent of if and only if v n is 

linearly independent of v^,...,v 

Using the first principle, we see that Ab = a b where a is 

m m n ,m n 


th 


and b is 

,, th 

the n 

n ,m 

i G 

1. Ab . = 

. E n a 

B 

l 

H 

O 

=0 n 


m 


■1 n-i r m-i 0 


For an 


. k, 
A b 


• = •£ « 
m-i 0 i]=0 


a . b 

n-i, n-i., m-i, 


. k, 
A b 


i 


m-i, 


Thus A b 


x Q n 

» ■ • • 


m+i. 


ll=°- ■ -i 2 =o 


1 

• • » Zj j-k 3. 

i,.=0 n-i, n-i, ,n-i. 


m-i. 


k v 11 k 1 k k-1 
is a linear combination of columns of A 


n-i» , n- 


• >— / 

h n 


-l n ,m-l , 


Using the second principle, Ab m is linearly independent of 

b ,...,b. if and only if a is linearly independent of b , ...,b. . 
ml n ml 



Thus for purposes of checking independence Ab^ can be represented 

by a n » Similarly A b m _^ can be represented by the leftmost column 

of A in the linear combination which has a nonzero coefficient. 

Zero coefficients are caused by previous vectors being dependent. 

Thus A b, . can be represented by column a r /i i \ , . , where 

m-i Q c J n-[(k-l)m + lgl+h 

h is the number of vectors already found to be dependent on previous 

k k-1 

vectors. To see this, consider A b . = A (A b . ) . The index 

m-i 0 m-1 0 

of the column representing A b will be the first (from the top) 

m 1^ 

k-1 u 

nonzero element of A b . . Since we move up a row for every 

m-i 0 

k-1 

independent vector, the first nonzero element of A b has index 

m i Q 

n-[(k-l)m + ig3 + h. 
k 

Thus A t> m _^ can be dependent on the previous vectors in two 

k-1 

ways. One way is that A b m _i was dependent on vectors previous 


to i.t. The second way is that the column a 


n- [ (k-1) m + i^ ] + h 


is 


dependent on columns to the right of it. In this second way 


A b„ . is the lowest power of A times b . which is dependent 

m-i 0 * 


on previous vectors. Thus there is a one to one correspondence 
between the ordered sets of Kronecker indices (ordered by def 3.3) 
and the zero patterns of the block triangular form. | | 

Once we have the block triangular matrix H, how do we retrieve 
the lead variables (to be defined) for the Brunovsky form? To see 
how this is done and to better understand block triangular systems, 
it is useful to think, in terms of derivative levels. 


Definition 3.6. Let x = Ax + Bu be a control system in block 



triangular form. A state variable x^ is said to be on the j ^ 
derivative level if it takes j derivatives of x^ to reach the 
controls. For instance in the block triangular system illustrated 
in equations (3.4) and (3.5), Xg,x^, and x^g are on the first derivative 
level. Xg,Xg, and x^ are on the second derivative level, x^ and 
are on the third level, and x 2 is on the fourth level. Lastly 
the variable x 1 is on the fifth derivative level. 

Definition 3.7 . Let x = Ax + Bu be a control system in block 
triangular form. Let x^ be a state variable on the derivative 

level. Then ^ is said to be a lead variable if it cannot be 
reached by taking the derivative of a state variable on the (j+l)th 
level. 

For all block triangular systems, x^ is a lead variable. 

For the particular system in equations 3.4, x 4 and x ? are also 
lead variables. Notice that x^ is on the third level; x 2 is the 
only variable on the fourth level but its derivative is not a 
function of x^ . Also x^ is on the second level; x^ and x^ are on 
the third level, but their derivatives Ulo not involve x^. 

Theorem 3.2. Let x = Ax + Bu be a control system in block 
triangular form. The derivative levels of the lead variables 
are precisely the Kronecker indices. 

Proof: This theorem is really just a restatement of Theorem 3.1. 

Notice that a lead variable occurs because of the inability to 
reach the variable from a higher level which in turn is caused 
by a column dependency in A. j | 

Once we have the lead variables in the block triangular form, 
we simply let them be the lead variables in the Brunovsky form. 
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For instance for the block triangular system x = Ax + Bu of equation 

y 6 ~ x 4 y 9 x 7 

y 7 - h y 10 = y 9 

y 8 = *7 v 3 “ *10 

v 2 - ?8 

We know the respective lengths of the integral strings because 
we know the derivative levels of x^,x 4 , and x^. 

Again, the essential characteristic of the block triangular 
system is that the controls do not appear "too soon". Thus the 
state variables in Brunovsky form are functions only of the state 
variables in the original system (that is, not functions of the 
controls) . 

Theorem 3.3. The above transformation from block triangular form 
to Brunovsky canonical form consists of only coordinate changes 
and linear feedback. 

Proof: As in the single control case, the derivative of x^ is a 

linear combination of x . 1 s except when x^ is on the first 

derivative level. In that case is a linear combination of 

the Xj's and u^'s. Subsequent derivatives of the state variables 

are also linear combinations of the appropriate variables. Let 

x i be a lead variable and let k be the derivative level of x. . For j < k, 

the j derivative of with respect to t,x^ ^ , is a linear 


3.4, we let 
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(k) 

combination of x 's and x. is a linear combination of the 

36 X 

x^'s and u^'s. Since the lengths of the strings in the Brunovsky 
form are determined by the derivative levels of the lead variables, 
our proof is complete. Q 

IV. Algorithms 


The actual algorithm for finding a matrix C to transform a 
matrix H to block triangular form is a generalization of the 
method found in Stewart [9] for transforming a matrix to Hessenberg 
form. 

We do this by placing the appropriate zeros in H one column 
at a time starting with the rightmost column. For each column, we 
multiply on the left and right by what is called an elementary 
reflector or Householder tranformation u, (see Stewart [9]). Each 

m K 

. 2vv x 

of the u^ is of the form u^ = I - -tj where I is the identity 

matrix and v is a vector. To illustrate this process, suppose we 
want to zero out the first k-1 elements in some column and at the 
same time insure that the k^ element is nonzero, suppose 
N is the column in question then we multiply 




0 


'*+ II v 
0 

where ||v|| is the Euclidean norm of 


k-1 


ll v ll ± vk 1 1 v II _ 
II V II (v + II V II ) 


0 

-II V || 


V, 


Starting with an index = m and k = n+m we zero out the first k-index 


-1 in elements . column k... . Then reduce k by 1 (or move one column to the 
left) and continue. If we find that the first k-index elements of 
column k are already zero, we reduce the index by 1, reduce k by 1 
and continue. All of the work is done by multiplication on the 
left, but it is easy to see that multiplication on the right does 


not undo the work. For instance u v , the matrix which produces 
th 


where C is at 


zeros in the k 1 "” column is of the form u. = S ^ 

k IP Ij 

most r (k-l) x (k- 1) . Thus multiplication of H by u^ on the right cannot 
affect columns k thru n+m, precisely the columns that have already 
been transformed. 

We give the essential part of the algorithm below. In so 
doing we use the pseudo code INFL of Stewart [4]. Here N is the 

dimension of the system and M is the number of controls. 


1. ) INDEX ■= M 

2. ) For K = N+M, ... ,1 

1. ) If H (K+l - INDEX , K +1) =0 

1) INDEX = INDEX -1' 

2. ) If K-index £ 1 Exit Loop 

3. ) ETA = max { |H(I,K) | : I = 1, . . . ,K- INDEX} 
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4. ) If ETA = 0 Step K 

5. ) V ( I , K) = H ( I , K) /ETA I = 1 K -INDEX 

6. ) SIGMA = sign (V(K - INDEX, K) ) / V 2 ( 1 , K) + . . . +V 2 (K-INDEX , K) 

7. ) V (K-INDEX, K) = V (K-INDEX, K) + SIGMA 

8. ) PI ( K ) = SIGMA * V (K-INDEX, K) 

After finding the vector V and the scalor tt. , the elementary 

T 

vv 

reflector U consists of U - I - - and the transformation 

k k w 

matrix C is C = u i U 2 ** ,u n + m * course some of the may be 

identities . 

A very similar procedure was used by Minimis and Paige [5] , 
as a first step, for the purpose of placing eigenvalues. They use 
a generalized upper Hessenberg matrix rather than a generalized, 
lower Hessenberg matrix. With renaming of coordinates they would 
have obtained the block triangular form. Once in block triangular 
form, as we have shown in this paper, the transformation to 
Brunovsky canonical form and hence the placing of eigen values is 
extremely easy. In Brunovsky form, the system is decoupled into 
several single control systems, each represented by a string of 
integrators . 

Once we have a block triangular system, we can utilize the 
zero pattern of H to select the lead variables for the Brunovsky 
form. We do this by working our way up the super diagonals of H 
until encountering zeros. The column number of the first zero 
in each of the first m superdiagonal working upwards tells us 
which variables have been 'kkipped over" and must be lead variables. 

Recall also that x^ is always a lead variable. The following 
algorithm shows precisely how this is done. 



1.) K = M 


2. ) I = N 

3. ) For L = 1, N +M 


1.) 

If 1 

H ( I , I +K) | fi 0 


1.) 

If I = 1 Exit Loop 


2.) 

I = 1-1 


3.) 

Step L 

2.) 

IOTA (K) = I+K 



N 

3.) 

Z(K) 

= Z C (I+K , J) X(J) 



J=1 

4.) 

K = 

K-l 

5.) 

Step 

' L 

If 

|H(1, J+K) | t 0 

1.) 

For 

t* 

II 

H 1 


1.) 

IOTA (L) = L 



N 


2.) 

Z (L) Z C (L, J) X(J) 


J=1 

"r Notice the lead variables Z of the Brunovsky form are computed 

in terms of the original x, not in terms of the intermediate block 

N 

triangular system. The variable Z (K) = Y(I+K) = Z C(I+K,J) X(J) 

J=1 

would be the I+Kth variable in the block triangular system. 

Once we know the lead variables in the Brunovsky form, all 
we need to know are the Kronecker indices. These can easily be 
found by counting the number of derivatives necessary to go from 
the lead variables in block triangular form to the controls. The 
following algorithm does this, using the value of IOTA from the 
previous algorithm. 



1.) For K=l, M 

1. ) I = IOTA (K) 

2. ) KAPPA (K) = 1 

3. ) For J = N+M, ... ,1 

1. ) If |H(I,J) | / 0 

1. ) If J > N Step K 

2. ) I = J 

3. ) KAPPA (K) = KAPPA' (K) + 1 

4. ) Begin J Loop again 

2. ) Step J 

4 . ) Step K 

V. RESULTS 


7 

We apply our theory and algorithms to a linear system on 3R 
with 3 controls. 

Example: The following results were achieved using single 

precision Fortran on a VAX 11/780 machine. For purposes of 
distinguishing zeros, numbers having absolute value less than 
0.00001 were considered zero. 


Original System 


*1 


o 

1 

0 

1 

0 

0 

0 


X 1 


0 

0 

I 


CM 

•X 


0 

0 

1 

2 

0 

0 

1 


X 2 


0 

0 

0 


X 3 


0 

0 

0 

1 

o ' 

0 

0 


x 3 


0 

6 

0 


~ U 1 

*4 

= 

0 

0 

0 

0 

1 

0 

2 


X 4 

+ 

0 

0 

0 


U 2 

*5 


0 

0 

0 

0 

0 

2 

0 


X 5 


0 

0 

0 


u 3 

*6 


0 

0 

0 

0 

0 

4 

5 


X 6 


0 

3 

0 


X 7 


3 

2 

0 

0 

0 

1 

6 

J 


X 7_ 


2 

1 

0 
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Transformation to Brunovsky form 

Z 1 = 0 . 16273 X '2 + 0.16781 x 3 + 0.75770 x 4 + 0.50853 x 5 - 0.33563 x g 
= 3 

Z 4 = - 0.85120 x 2 - 0.09259 x 3 - 0.03736 x 4 + 0.48083 x g + 0.18518 x g 


Here k^,ic 2 , and k 3 are in order as in definition 3.1 and are not 
necessarily the Kronecker indices associated with the original 
controls u^u^ and u 3 respectively. That is, the subscript 1 does 
not identify with u^. In a Brunovsky canonical form, the 
controls appear only at the end of integral strings; that is , after 
three derivatives of Z^, three derivatives of Z 4 , and one derivative 
of Z^. All state variables Z should be functions only of the state 
variables x. 

To get some idea of the accuracy of our numerical method we 
will look at the actual results in the first string. 

We have 

m 

z 2 = Zf = 0.16273 x 3 + 0.49327 x 4 0.75770 x 5 - 0.32546 x g 
- 0. '00002 x^ - 0.00003 u 2 

Notice that the control u 2 does appear in Z 2 , but only negligibly. 
Computing, we find 

Z_ = Z„ = Z* = - 0.00060 x. - 0.00004 x„ + 0.16273 x. + 0.49327 x c 
3 2 1 1 2 4 5 

+ 0.01354 x c - 0.64098 x., - 0.00004 u, - 0.00002 u- 

6 7 12 

- 0.00003 u 2 



We again have in Z. that the appearance of the controls u, and u 

J 12 

and the derivative u 2 is negligible. 

* r* - • » « 

The new control = z 3 = z 2 = Z i = 

- 1.92294 x x - 1.28202 x 2 -0.00004 x 3 - 0.00014 x 4 

+ 0.16273 x 5 + 0.39972 x g - 3.45276 x ? 

- 1.28196 - 0.60036 u 2 - 0.00024 u 3 - 0.00004 fc 

- 0.00002 u 2 - 0.00003 u 2 

In the new control v. , the original controls appear substantially 
but derivatives of controls appear only negligibly. These results 
can be improved by using more precise arithmetic, but we want to 
show how well it works even in single precision. 

The weakest link in the algorithm is of course the recognition 
of zeros. Because the computer uses discrete arithmetic, there 
are no absolute zeros and nonzeros. We must have some criteria 
for deciding which numbers are zero. We chose ah arbitrary cutoff 
point of 0.00001 in our example. Perhaps a more objective approach 
would be to compare an element with other elements in the matrix 
by using a matrix norm. We expect however that in an actual 
application, there will always be some "tailoring" of this cutoff 
point to suit the particular situation. 

As mentioned previously, George Meyer is using transformations 
to Brunovsky canonical forms to design automatic flight controllers 
for vertical and short take off aircraft. In his scheme, the weak 
dependencies of the new state variables on the controls in our 
example are ignored, and errors are treated by a regulator which 
resides -in the aircraft's on board computer. 



In practice, we may be receiving noisy measurements for 
state and control variables. For our process, once we reach block 
triangular form, time derivatives are required to complete the 
transformations to Brunovsky form. Meyer [2] has a beautiful 
technique for moving from block triangular to Brunovsky form which 
involves "smoothing integration" using the inflight computer before 
the dif ferentation process takes place. 

The theory of transformations of nonlinear systems to linear 
systems is developed in [11] , [12] , [13] , [14] , [15] , [16]. Recent 

applications of this theory to automatic flight control are found 
in [17] , [18] , [19] . The techniques of this paper are employed in [8] 
to build approximate transformations of nonlinear systems to 
Brunovsky form. 
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